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Abstract 



We calculate the cross sections for the production of one and more electron- 
positron pairs due to the strong electromagnetic fields in relativistic heavy 
ion collisions. Using the generating functional of fermions in an external 
field we derive the iV-pair amplitude. Neglecting the antisymmetrisation in 
the final state we find that the total probability to produce N pairs is a 
Poisson distribution. We calculate total cross sections for the production 
of one pair in lowest order and also include higher-order corrections from 
the Poisson distribution up to third order. Furthermore we calculate cross 
sections for the production of up to five pairs including corrections from the 
Poisson distribution. 
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I. INTRODUCTION 



Currently new accelerators are being built at CERN ("Large Hadron Collider", LHC) 
and at Brookhaven ( "Relativistic Heavy Ion Collider", RHIC), which allow heavy ions up to 
lead and uranium to be accelerated to highly relativistic energies. The main motivation for 
this is to advance into regions of baryonic densities and temperatures in order to observe the 
formation of the quark gluon plasma. But this has also renewed the interest in studying pure 
electromagnetic processes in these collisions, especially the electron-positron-pair production 
in peripheral collisions, where both ions do not interact with each other. 

On the one hand pair production in central collisions is a signal one wants to use in 
order to study the properties of the quark gluon plasma [TJ,|2]]. The pure electromagnetic 
pairs are a large background process to these. On the other hand these processes can also 
be studied already in existing accelerators ("Super Proton Synchroton", SPS) M. Recently 
anti-hydrogen atoms were produced at LEAR ( "Low Energy Antiproton Ring" ) in a similar 
process, where antiprotons in collisions with Xenon produce electron-positron pairs, where 
the positron is not produced as a free particle, but is captured in a bound state of the 
antiproton Q. 

The theoretical treatment of this problem goes back to the beginning of QED 0. But 
it was remarked only recently that calculations in lowest order perturbation theory violate 
unitarity |J . In the meantime several authors have studied this problem and found that this 
means that the production of multiple pairs becomes important then. All found that the N- 
pair creation probability can be written approximately in the form of a Poisson distribution, 
which then restores unitarity [0-|9|]. It was shown that generally the iV-pair amplitude 
can be reduced to the product of the vacuum amplitude and the antisymmetrised product 
of N reduced one-pair amplitudes. The neglect of exchange terms in the calculation of the 
probability then leads to a Poisson distribution |10| . Corrections to this Poisson distribution 
were calculated in the case of two-pair production and were found to be small (around 1%), 
its importance decreasing also with higher Lorentz factors. 

Whereas large impact parameter contribute to the single-pair production, the iV-pair 
production with increasing iV is dominated mainly by small impact parameters, smaller 
than the Compton wavelength of the electron. Thus the "equivalent-photon-approximation" 
(EPA) cannot be used anymore but exact calculations of the reduced probability are needed, 
since the momentum of the virtual photon is not small compared to the electron mass. 

In Sec. II we show how the iV-pair production amplitude can be found using the path 
integral formalism. We find again the results of ||10|| . We then discuss in Sec. Ill our 
approach to the calculation of the cross sections for one- and multiple-pair production taken 
into account only peripheral collisions. Results for this are given in Sec. IV and compared 
with EPA. Differential cross section for the one-pair production are shown too. 



II. iV-PAIR PRODUCTION IN AN EXTERNAL FIELD 



In the path integral formalism the generating functional for a fermion field ip in a external 



field can be written as |11 



Z[rj, fj,A) — N / Vi/jVi/j expli / d 4 x 



C(x) 
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with the Lagrangian 



+r)(x)ip(x) + i/j(x)rj(x) 



C(x) = ijj(x)(i Jp — m)il>(x) 



(1) 



(2) 



where = + ieA^ is the usual gauge invariant derivative. We normalize Z[rf,fj,A] so 
that it corresponds to the identity in the absence of the fermionic sources and external field. 
We separate the fields tp = if) c i + if)' and if) = ifjci + tfj' into their classical parts and fluctuations 
around them. Using Feynman boundary conditions in order to invert the equations of motion 
the classical solutions are given by 



d 4 y S A (x,y)rj(y) 



ipd(x) = ~ I d 4 yri(y)SA(y,x), 
where Sa is the Dirac propagator in the external field 

(i p-eJL- m)S A (x, y) = 5 4 (x - y). 
With this we can rewrite the functional in Eq. (JJ) in the form 

det - m)\ 



(3) 



(4) 



det — m)) 



Z[r},ri,A] 



(5) 

The Green functions of the Dirac fields ip and ^ can now be found by varying Z with respect 
to the sources 



i x n) y-a) 
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5ij(y n ) 5r](x n 



-Z[rj,rj,A] 
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We use the asymptotic out-system. The initial and final states for iV-pair creation are then 



\z, out) = |0) 

\f, out) =dt utsi (p 1 )b+ uts ,(p' 1 ).. 



d 



out sjv 



(pN)K uts ,(p' N )\o) 



(7) 



where the one-particle creation operators b + and d + obey the usual anticommutation rela- 
tions for free particles. In order to get the amplitude for the creation we make use of the 
Lehmann-Symanzik-Zimmermann reduction |12[|. The first term in Eq. fl5|) can be identified 
with the vacuum amplitude [|K|. Thus we get the amplitude in the form 
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(f\S\i) = (0\S\0) J d i x 1 . . . d 4 x N d* yi . . . d^y N 

N _^ 
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(8) 



The variations in Eq. (H) lead only to connected graphs. We extract the individual fermion 
lines with the help of the Wick theorem [Tj| in order to get 



(f\S\i) = (Q\S\0) £ e ff 5V»^(i))-'5V.M«)) 



(9) 



where we have summarized momentum and spin quantum numbers of electrons as Ui = 
and of positrons as Vj = (pj,Sj) and summed over all permutations of Sn- The 
reduced one-pair amplitude S R (ui,Vj) describes a single fermion line interacting with the 
external field to arbitrary order. They are given by expanding the Dirac Propagator Sa and 
doing the integrals in Eq. (||), see also |10| , 



n=2 



d i z 1 ...d i z n u { <\p' i )e ip '^ 
x|[e 4(z n )) S F (z n - Zn-x) . . . [e 4{zx)) 



(10) 



The A^-pair amplitude in Eq. @ is therefore the product of the vacuum amplitude and an 
antisymmetrised product of the reduced amplitudes. Due to the correlation between the 
produced electrons and positrons it is justified to neglect the Pauli principle in the final 
state, i.e., terms in the probability for TV pairs coming from two different permutations are 
not taken into account. This leads finally to the probability to produce iV pairs 



P(N) 



N 



(11) 



where the vacuum amplitude shows up as the exponential function. As shown in the 
neglect of the exchange terms does not give rise to large deviations for total two-pair creation 
probabilities. Therefore the use of the Poisson distribution seems to be justified. With 
the knowledge of reduced one-pair probabilities, Eq. QTTJ) allows the calculation of A^-pair 
probabilities and A^-pair cross sections. 



III. CALCULATION OF CROSS SECTIONS 

The first order that contributes to the reduced one-pair amplitude is the second order. For 
the creation of an electron with momentum p_ and spin s_ and a positron with momentum 
p + and spin s + we get the amplitude 
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(2vr) 4 ' ^ L ' p 2 -m 2 
x4(p ++ p)\v^\p + ). (12) 



We describe the electromagnetic potential of the two ions as a superposition of the four- 
potentials 

F(q) 



A,(q) = -2nZe^{uVe^5( qU V) 



+ufe-^5(qu^)Y (13) 

using the straight-line approximation for the trajectories of the ions, b is the impact param- 
eter, w 1 ' = 7(1,0,0,/?) = jw^ 1 ' and = 7(1,0,0,—/?) = are the four-velocities of 
the ions in the center-of-velocity frame. Throughout the calculations we are making use of 
the monopole form factor 

F M = (14) 

with A = 83 MeV for Au ions, in order to describe the extended charge distribution of 
the ion, as this form factor can be treated analytically [T5 ]. One-pair cross sections are 



not sensitive to the detailed choice of the form factor. The cross sections for the multiple- 
pair productions are more sensitive to this as they are produced mainly at smaller impact 
parameters. 

We get the total reduced probability by summing the absolute value squared of Eq. ([12]) 
over all electron and positron spins and integrating over all momenta 

PH^ffffi-g-xW. (18) 

By rewriting the spin summation as a Dirac trace over 7-matrices we get 

1 
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(16) 



where tr is the usual trace over the Dirac indices. The 5 functions in the amplitude in 
Eq. (|l^) determine the longitudinal components of the integration variables. We still have 
to integrate over the components transverse to the velocity of the ions. We have changed 
these variables to the momentum of one of the photons q\ = p~ — p and to the difference 
q = q[ — q\. This allows us to rewrite the 6-dependent probability in Eq. fllED in the form 
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P R (b) = I d 2 qP H (q)e 



R( 



7 iqL 
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where the Fourier transformed probability is 
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with the propagators 
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The momenta ki and parameters of the propagators in Eq. 
0) 



(19) 

are defined as (see also 
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(20) 



The trace over the Dirac 7-matrices is done with the algebra-program FORM [16|. The 
advantage of calculating the Fourier transform of the total probabilities is due to the absence 
of oscillating terms in Eq. (|P7|), which would otherwise prevent convergence in the Monte 
Carlo integrations. 

From this we can get the total cross sections a by integrating the probabilities (Eq. ([TT|)) 
over the impact parameter 



a(N) = 2ir / dbbP(N,b), 

J2R 



(21) 



starting with 2R (for symmetric collisions) in order to take only peripheral collisions into 
account. 

We can use P(N = 1,6) directly as given in Eq. (|TTD. But for large impact parameter 
the integrand oscillates very fast allowing not a very accurate calculation of P. Therefore we 
choose to do the calculation for the cross section for the one-pair production by expanding 
the Poisson distribution of a(N = 1) up to third order 



a(N = 1) = - 



3! 
2! 



a 



(3) _ 



(22) 
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Integrating from b = on the can be written as 

a% = J d 2 bd 2 qe^P R (q) 



= {2n f J d 2 q5 2 (q)P R (q) 

= (2n) 2 P R (0), (23a) 
4% = ^fjdqq{P R (q)} 2 , (23b) 

4% = ^fj P R (Qi)P R (<l2)P R (\qi + (23c) 

Therefore we can calculate the cross section by using the Fourier transformed P R (q) alone. 
This leads to more accurate results for the cross section as the cross section for the one-pair 
production is dominated by large impact parameter, where the impact parameter depen- 
dent probability is not very accurate due to the large number of oscillations in the Fourier 
transform. 

For the calculation of we have to set q to zero in Eq. (|i"8|) , i.e. q\ ' = q\. This leads 
to a number of simplifications 

N 3 = N , N 4 = N 1} N bD = N 2D , N 5X = N 2X . (24) 

The integrations over q\ can then be done analytically as they are standard two- 
dimensional Feynman integrals. The q\ 's in the numerator only appear in connection with 
the scalar products p+qi, p~q\ and q\qi and the fej's in the propagators are also either 
—p-, —p+ or their sum. Therefore by using the relations 

<fi = N 2 - ml 

2%.% = Ni- N -r t , (25) 

we can express the scalar products in the numerator as propagators and scalars in order to 
eliminate all qi's from the numerator. In the end we are left with scalar Feynman integrals 
with different numbers of propagator terms in the denominator. The higher integrals with 
quadratic propagator terms cannot be reduced to simpler integrals with the method de- 
scribed in ]n|,|T7| • But we solve them by differentiating the standard integrals with respect 



to the parameters m\ and m|. 

The integration over the momenta of electron and positron in the final state are done 



with the Monte Carlo integration routine VEGAS |L8]|. We write the momenta in polar 
coordinates in the plane transverse to the velocity of the ions. The integration over one 
of the angles is trivial and there remains a five dimensional integral over momenta and 
the relative angle between the transverse electron and positron momenta. The integration 
boundaries were incremented until the Monte Carlo integration converged with accuracies 
better than 0.5%. We did not take any form factors into account for P i? (0)'s as the results 



of [14 1 show that the form factor is only important for small 6's, i.e., only for large q's. At 



small q's differences between calculations with and without form factors are not observable. 
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For the second order correction (Eq. ( |23b|) ) we calculate the values of P R (q) for 
q 7^ with the method described in |l4yT5|l . We are adding to this the value at q — 0, and 
approximate the P R (q) by cubic splines. We use a logarithmic scale for q for the interpolation 
in order to get a better approximation for the steep increase at small g's. 

The calculation of the third order correction (Eq. ( |23c| )) follows that of a^. 

Finally we have to correct for the fact, that we are integrating from 6 = instead of 2R. 
First we calculate the impact parameter dependent probability (Eq. fllTD) by integrating 
over the angular part and get 

P R (b) = 2nJ dqqP R (q)J (qb), (26) 

with the Bessel function J Q . Then we subtract from the the contribution from < b < 
2Rby integrating 

*£. = 2.6™, (27) 
db n\ 

from b = to 2R. As P R (b) only decreases slowly with b, large impact parameters are 
still important for the cross section. The Bessel function in Eq. ( p6|) oscillates fast for 



large b, which would lead to a poor accuracy for if one would integrate over the impact 
parameter directly. 

For increasing N pair production occurs within increasingly smaller 6's as the external 
field has to supply the minimal energy of 2Nmc 2 to create N pairs, which further restricts 
the system. In this case the influence of the large 6's decreases and we can calculate the total 
multiple-pair cross sections exactly by using the whole Poisson distribution and integrating 
directly over the impact parameters 

a(N) = 2n J™ dbb [ ^f- Y p ». (28) 



IV. RESULTS AND CONCLUSIONS 

Using the method described in the last section we have calculated the impact parameter 
dependence of the one- and multiple-pair cross section up to three pairs. In Fig. [l| we show 
the results of the calculation with and without a monopole form factor for a center-of-mass 
Lorentz factor 7 = 100 and for a 197 Au 79 — 197 Au 79 collision. Also shown are the results of 
the equivalent photon approximation (EPA) using for the reduced one-pair probability in 
the Poisson distribution 

pEPA ^ = §i(Z*f(fi^ 2 (^\ (29) 

with 5 = 0.681. A c is the Compton wavelength of the electron and jtar = 2j 2 — 1 is the 
Lorentz factor in the target system. The EPA approximation of Eq. (BT3) is only valid for 



impact parameters A c 7f ar 5 > 2b > 2X C . We find in agreement with earlier results |14|T5 
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that the EPA gives results which are much to large at small b. As already discussed before, 
the one-pair production cross section decreases only slowly with larger impact parameter in 
contrast to the multiple-pair cross sections. 

In Fig. H we show the dependence of the total cross section for the one-, two- and three- 
pair production for a 208 Pb 82 — 208 Pb 82 collision on the Lorentz factor 7. We plot the results 
for the calculation with a monopole form factor using the full Poisson distribution and also 
using the lowest order alone. One sees that the Poisson distribution reduces the multiple- 
pair cross sections considerably whereas it seems to have only a small effect on the one-pair 
production. Also shown are the results of the EPA approximation by integrating only over 
impact parameters larger than the Compton wavelength. They are in fair agreement with 
the exact results, even though they tend to be either to large or to small systematically. 

Our results for the total cross sections are summarized for the different colliders in Tables 
I and [II[ The differences between the calculations with and without a monopole form factor 
show the dependence on the detailed form of the charge distribution in the ion. More realistic 
form factors are needed in order to describe pair production with a larger N more accurately. 
Moreover our results are in agreement with the results of fim , which were calculated only 
at lower 7 and with a Monte Carlo integration alone. 

Recently anti-hydrogen atoms were produced in a LEAR experiment at CERN by col- 
liding antiprotons with Xenon |3J with 7 tor = 1.8. With our approach we get a total cross 
section for the creation of a free electron-positron pair in this case of = Z 2 o~o, with 
o"o = 14 nb. This is more than three orders of magnitude larger than the bound-free transi- 
tion, which leads to anti-hydrogen formation. 

We have applied our calculations only to electron-positron-pair production. Nevertheless 
the pair production of heavy leptons fi and r gets more important with increasing energies of 
the ions. With our method we can calculate cross sections for this processes, too, considering 
only peripheral collisions. Moreover only single-pair cross sections are important here due 
to the larger masses. 

Finally we also give differential cross sections for the energy and the angle with the beam 
axis of either electron or positron. We get these by integrating the first order of the one-pair 
cross sections in Eq. (^) with the help of the Monte Carlo integration routine sorting the 
results into bins. The results are shown in Figs. |3] and |j. Most electrons are produced rather 
close to the beam axis and also with relatively low energy of the order of a few MeV. 
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FIG. 1. The b dependent differential cross sections (ft) for the iV-pair production in a 

19 Au 79 — 197 ^4-u 79 collision with 7 = 100 for iV = 1,2,3. The solid lines show the results for the 
calculation with a monopole form factor, the dotted lines for a point charge and the dashed lines 
are the EPA results. 
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FIG. 2. The 7 dependent total cross sections <r(iV)(7) for the iV-pair production in a 
208p^82 _208 p^82 collisions, for N = 1,2,3. Exact calculations are shown as solid lines, dashed 
lines are the results for the corresponding first orders 0-^(7) and the dotted lines are the EPA 
results. 
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FIG. 3. Differential cross sections 
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FIG. 4. Differential cross sections 
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TABLES 







with 


monopole form factor (kb) 






7 


a(N = 1) 


a(N = 2) 


a(N = 3) 


a(N = 4) 


a(N = 


5) 


8.916 (SPS) 


2.29 


9.89 10~ 2 


6.26 10" 3 


4.52 10" 4 


3.2510 


-5 


100 (RHIC) 


31.8 


0.899 


0.1122 


1.827 10" 2 


3.1110 


-3 


3400 (LHC) 


223 


4.00 


0.785 


0.219 


6.9410 


-2 



TABLE I. Results of the calculation for the production of N = 1, . . . , 5 pairs with a monopole 
form factor. For SPS and LHC we show total cross sections for 208 Pb 82 — 208 P6 82 collisions and for 
RHIC total cross sections for 197 Au 79 — 197 ^4u 79 collisions. 







without monopole form fact 


or (kb) 
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a(N = 1) 


a(N = 2) 


a(N = 3) 


a{N = 4) 


a(N = 


5) 


8.916 (SPS) 


2.28 


0.1019 


6.66 10" 3 


5.08 10" 4 


3.9510 


-5 


100 (RHIC) 


31.8 


0.897 


0.1136 


1.881 10" 2 


3.3710 


-3 


3400 (LHC) 


223 


3.95 


0.776 


0.219 


7.0910 


-2 



TABLE II. Results for the calculations for the production of = 1, ... ,5 pairs for a point 
charge. The ions used for the different 7's are the same as in Tab. §. 
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